Massive, massless and ghost modes of gravitational waves from higher-order gravity 
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We linearize the field equations for higher order theories that contain scalar invariants other than 
the Ricci scalar. We find that besides a massless spin-2 field (the standard graviton), the theory 
contains also spin-0 and spin-2 massive modes with the latter being, in general, ghost modes. Then, 
we investigate the possible detectability of such additional polarization modes of a stochastic gravi- 
tational wave by ground-based and space interferometric detectors. Finally, we extend the formalism 
of the cross-correlation analysis, including the additional polarization modes, and calculate the de- 
tectable energy density of the spectrum for a stochastic background of the relic gravity waves that 
corresponds to our model. For the situation considered here, we find that these massive modes are 
certainly of interest for direct detection by the LISA experiment. 
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I. INTRODUCTION 

Recently, the data analysis of interferometric gravita- 
tional wave (GW) detectors has been started (for the 
current status of GWs interferometers see pHII) and the 
scientific community aims at a first direct detection of 
GWs in next years. The design and the construction of a 
number of sensitive detectors for GWs is underway today. 
There are some laser interferometers like the VIRGO 
detector, built in Cascina, near Pisa, Italy, by a joint 
Italian-French collaboration, the GEO 600 detector built 
in Hannover, Germany, by a joint Anglo-German collabo- 
ration, the two LIGO detectors built in the United States 
(one in Hanford, Washington and the other in Livingston, 
Louisiana) by a joint Caltech-MIT collaboration, and the 
TAMA 300 detector, in Tokyo, Japan. 

Many detectors are currently in operation too, and 
several interferometers are in a phase of planning and 
proposal stages (for the current status of gravitational 
waves experiments see 6-8]). The results of these de- 
tectors will have a fundamental impact on astrophysics 
and gravitational physics and will be important for a bet- 
ter knowledge of the Universe and either to confirm or 
rule out the physical consistency of General Relativity or 
any other theory of gravitation [9(. Several issues com- 
ing from Cosmology and Quantum Field Theory suggest 
to extend the Einstein General Relativity (GR), in order 
to cure several shortcomings emerging from astrophysi- 
cal observations and fundamental physics. For example, 
problems in early time cosmology led to the conclusion 
that the Standard Cosmological Model could be inade- 
quate to describe the Universe at extreme regimes. In 
fact, GR does not work at the fundamental level, when 
one wants to achieve a full quantum description of space- 
time (and then of gravity). 

Given these facts and the lack of a final self-consistent 



Quantum Gravity Theory, alternative theories of grav- 
ity have been pursued as part of a semi-classical scheme 
where GR and its positive results should be recovered. 
The approach of Extended Theories of Gravity (ETGs) 
based on corrections and enlargements of the Einstein 
scheme, have become a sort of paradigm in the study 
of the gravitational interaction. Beside fundamental 
physics motivations, these theories have received a lot of 
interest in cosmology since they "naturally" exhibit in- 
flationary behavior which can overcome the shortcomings 
of standard cosmology. The related cosmological models 
seem realistic and capable of coping with observations. 
ETGs are starting to play an interesting role to describe 
today's observed Universe. In fact, the good quality data 
of last decade has made it possible to shed new light on 
the effective picture of the Universe. 

From an astrophysical point of view, ETGs do not re- 
quire finding candidates for dark energy and dark matter 
at the fundamental level; the approach starts from taking 
into account only the "observed" ingredients (i.e. grav- 
ity, radiation and baryonic matter); it is in full agree- 
ment with the early spirit of a GR that could not act 
in the same way at all scales. For example, it is possi- 
ble to show that several scalar-tensor and /(i?)-models 
(where / is a generic function of the Ricci scalar R) agree 
with observed cosmology, extragalactic and galactic ob- 
servations and Solar System tests, and give rise to new 
effects capable of explaining the observed acceleration of 
the cosmic fluid and the missing matter effect of self- 
gravitating structures without considering dark energy 
and dark matter. For comprehensive reviews on the ar- 
gument, see [Toj . 

At a fundamental level, detecting new gravitational 
modes could be a sort of experimentum crucis in order 
to discriminate among theories since this fact would be 
the "signature" that GR should be enlarged or modified 
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The outline of the paper is as follows. In Sect. [H] the 
general action of the class of theories under considera- 
tion is introduced. Then we will linearize them around 
a Minkowski background to find the modes of the met- 
ric perturbations. In Sect. IIII1 we take into account the 
various polarizations of the massless and massive modes, 
while in Sect. IIVI we investigate the response of a sin- 
gle detector to a GW propagating in certain direction 
with each polarization mode. In Sect. [V] we discuss the 
spectrum of the GW stochastic background where also 
further modes are considered. Conclusions are drawn in 
Sect. EH 



II. HIGHER ORDER GRAVITY 



Expanding the third term on the RHS of (|2.5[) and 
using the purely geometrical identity G ab . b = 0we get: 

[2/ -RF- 2R ab V a V b (f P + 2/ Q ) - RD(f P + 2f Q ) 
-2(f P P + f Q Q)} (2.6) 

If we define 



$ = F+-(f P + f Q )R 
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and 

= RHS of (EES 



(2.7) 



Let us generalize the action of GR by adding curvature 
invariants other than the Ricci scalar. Specifically, we 
will consider the action 1 



5* = 



where 



J d 4 xV^f(R,P,Q) 

P = R ab R ab 

r\ — r> T>abcd 



(2-1) 



(2.2) 



Varying with respect to the metric one gets the field 
equations jl3l |: 



then we get a Klein-Gordon equation for the scalar field 



(2- 



In order to find the various modes of the gravity waves 
of this theory we need to linearize gravity around a 
Minkowski background: 



$ = $ + 



Then from eq. (|2.7|) we get 



(2.9) 



FG^ = ~g flv (f-RF)-(g fiU n-V fi V„)F 

abc \ 



'2 (fpR^Rau + fQ RabctiR ahC v ) 

-g^V a Vb(fpR ab ) - a(fpR^) 
-2V a V b (/ P + 2/ Q R a { ^ 



F : 



ff 



d£ , = df_ 
dP 1 lQ dQ 



(2-3) 
(2.4) 



where we have set 
df_ 
dR 1 

and □ = g ab V a^b is the d'Alembert operator while the 
notation T^j) = ^(Tij+Tji) denotes symmetrization with 
respect to the indices 

Taking the trace of eq. (|2.3p we find: 

. f F + y-i? 



i (2/ -RF- 2V a V b ((/ P + 2f Q )R ab ) - 2(f P P + f Q Q)) 

(2.5) 



S^ = dF+ 1(8 f P + 6f Q )R + l(f PQ + f Q0 )SR (2.10) 

where Rq = R^v) = and similarly f P0 = §p\, hlv 
(note that the indicates evaluation with the Minkowski 
metric) which is either constant or zero. By 5R we denote 
the first order perturbation on the Ricci scalar which, 
along with the perturbed parts of the Riemann and Ricci 
tensors, are given by (see for example Ref . [Tij ) : 



SR^u = - {dcrduh' 7 ^ + dadfjfv - dpd u h - Dh^) 
SR. = d^duh^ -Uh 
where h = rj^h^. The first term of eq. (12.101) is 
dF OF dF 

SF = lm^ R+ 8P^ p+ 8Q^ Q (2 - u) 

However, since SP and 5Q are second order we get 5F ~ 
F R q 6R and 
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S<S> = \F R0 + -(fp + f Q0 ) \ SR (2.12) 

Finally, from eq. (|2.6p we get the Klein- Gordon equation 
for the scalar perturbation (5$ 
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as<f> = - 



3 F tR0 + 3 (/po + /go) 
|jii a6 a a a b (/po + 2/ Q0 ) - yBD(fpo + 2/qo) 



(2.13) 



Using equations (I2.16p ~ d2.18p into (|2.15l) and after some 
algebra we get: 



fiv 



-ffc 2 -fc 4 /p0 t 4/Q0 ^ 
2 V -^o 



(2.19) 



The last two terms in the first line are actually are zero 
since the terms fpo, /qo are constants and we have de- 
fined the scalar mass as m 2 = — 



" s ~ 3 F, a „ + T(f PO +fQo) ' 

Perturbing the field equations (|2.3p we get: 
F {8R^ - -rin V 6R) = 

-(n^n - d»d v ){8<h - |(/p + f Q0 )SR) 
-V^dad b (f P0 SR ab ) - nifpoSR^) 
+2d a d b (f P0 5R\/ V) + 2/qo M2V) 6 ) 

(2.14) 

It is convenient to work in Fourier space so that for ex- 
ample d^h^y — > ik^h^y and Oh^ — > —k 2 h^ v . Then the 
above equation becomes 

FoiSR^ - ^p^SR) = 

(V^k 2 - k^k u ){8^ - ~(/ F0 + f Q0 )5R) 

+V^k a k b (f P0 SR ab ) + k 2 {f P0 SR^) 
-2k a k b (f P0 SR a (fl S%) - Ak a k b (f Q0 8R\^) 



(2.15) 



We can rewrite the metric perturbation as 



hfw = h^ v - - r)^ + rj^hf (2.16) 

and use our gauge freedom to define to demand that the 
usual conditions hold d^h^ = and h — 0. The first of 
these conditions implies that k f Ji llv — while the second 
that 



h = Ah f 



(2.17) 



Defining hf 
bations: 



we find the equation for the pertur- 



k 2 + 



spin2 



h^u = 



(2.20) 



where we have defined m^„,„ 9 = — — , while from 

spint Jpo+4/qo 



eq. ([2~T3"jl we get: 



Uh 



f 



(2.21) 



From equation (|2.20p it is easy to see that we have a mod- 
ified dispersion relation which corresponds to a massless 
spin-2 field (k 2 = 0) and a massive spin-2 ghost mode 



k 2 = t- 



F 



with mass m 2 



To see 



this, note that the propagator for h^ v can be rewritten 
as 



G(k) oc 



1 



k 2 k 2 



(2.22) 



spin2 



Clearly the second term has the opposite sign, which in- 
dicates the presence of a ghost, and this agrees with the 
results found in the literature for this class of theories 

Also, as a sanity check, we can see that for the Gauss- 
Bonnet term Cgb = Q — 4P + R 2 we have fpo — —4 and 
/qo = 1. Then, equation (I2.20p simplifies to fc 2 /j M „ = 
and in this case we have no ghosts as expected. 

The solution to eqs. (|2.20p and (|2.2ip can be written 
in terms of plane waves 

hpv = A^(~f) ■ exp(ik a x a ) + cc (2.23) 



hf = a(p) ■ exp(iq a x a ) + cc 



(2.24) 



where 



With these in mind we have: 



1 



SRpu = ^ (ZkfiKhf + k 2 ri^h f + k 2 h^ v ) 
SR = 3k 2 h f 

k a k/3 SR a ^ = -- ((k 4 !]^ - k 2 k^k v )hf + fc 4 /^„) 



k a k b 5R a ,J b u) = -k 2 k^Kh f 



(2.18) 



k a = {oj mBpm2 ,l?) 



z 2 
spin 



J a = (ui ms , p) 0J ms = ^m 2 s +p 2 



(2.25) 



and where m sp i n 2 is zero (non-zero) in the case of mass- 
less (massive) spin-2 mode and the polarization tensors 
^Vy(f^) can be found in Ref. [l8[ (see equations (21)- 
(23)). In eqs. (ET2"0)) and (ET2"3")l the equation and the so- 
lution for the standard waves of General Relativity [26[ 
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have been obtained, while eqs. (|2 . 2 1 1) and (|2.24l) are re- 
spectively the equation and the solution for the massive 
mode (see also |27|). 

The fact that the dispersion law for the modes of the 
massive field hf is not linear has to be emphasized. The 
velocity of every "ordinary" (i.e. which arises from Gen- 
eral Relativity) mode is the light speed c, but the 
dispersion law (the second of eq. (|2.25[) ) for the modes of 
hf is that of a massive field which can be discussed like 
a wave-packet [27J. Also, the group- velocity of a wave- 
packet of hf centered in p is 



vg 



1 



(2.26) 



which is exactly the velocity of a massive particle with 
mass m and momentum p. 

From the second of eqs. (|2.25l) and eq. (|2.26l) it is 
simple to obtain: 



III. POLARIZATION STATES OF 
GRAVITATIONAL WAVES 

Considering the above equations, we can note that 
there are two conditions for eq. (|2.13[) that depend on 
the value of k 2 . In fact we can have a k 2 = mode that 
corresponds to a massless spin-2 field with two indepen- 
dent polarizations plus a scalar mode, while if we have 
k 2 ^ we have a massive spin-2 ghost mode and there 
are five independent polarization tensors plus a scalar 
mode. First, lets consider the case where the spin-2 field 
is massless. 

Taking p in the z direction, a gauge in which only 
A\\, ^22, and A12 — A 2 \ are different to zero can be 
chosen. The condition h = gives An = —A 2 2- In this 
frame we may take the bases of polarizations defined in 
this way 2 



1 
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VG 



\J ij 2 — m 2 



(2.27) 



Then, wanting a constant speed of the wave-packet, it 
has to be US 




(3-1) 



Now, putting these equations in eq. (|2.16p . it results 



m = V i 1 ~ v gP 



(2.28) 



Now, before we proceed with the analysis, we should 
discuss t he p henomenological limitations to the mass of 
the GW [Hf. Taking into account the fact that the GW 
needs a frequency which falls in the range for both of 
space based and earth based gravitational antennas, that 
is the interval 1(T 4 Hz < f < lOKHz [U-d Hj|3l|, a 
quite strong limitation will arise. For a massive GW, 
from [32| it is: 



2nf = uj = ^m 2 +p 2 , (2.29) 
were p is the momentum. Thus, it needs 



OeV < m < W- n eV. 



(2.30) 



A stronger limitation is given by requirements of cos- 
mology and Solar System tests on extended theories of 
gravity. In this case it is 



h^{t,z) 



A+(t-z)e^ + AX(t-z)e^ 



h a (t 



v G z)e s 



(3.2) 



The terms A+(t 



z)e {+ 



+ A x {t - z)e 



(x) 



describe 



the two standard polarizations of gravitational waves 
which arise from General Relativity, while the term 
h^t—VGZjri^u is the massive field arising from the generic 
high order f(R) theory. 

When the spin-2 field is massive, we have that the bases 
of the six polarizations are defined by 



=,(+) 



1 



1 
..0-10 
v 2 \ 



1 



1 



7! 

V 2 \ 1 



(x) 



,(C) 




»W = 




OeV < m < \0- 33 eV. (2.31) 

For these light scalars, their effect can be still discussed 
as a coherent GW. 



2 The polarizations are defined in our 3-space, not in a spacetime 
with extra dimensions. Each polarization mode is orthogonal to 
one another and is normalized e M „e'"' = 2<5. Note that other 
modes are not traceless, in contrast to the ordinary plus and 
cross polarization modes in GR. 
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Figure 1: The six polarization modes of gravitational waves. 
The picture shows the displacement that each mode induces 
on a sphere of test particles at the moments of different phases 
by 7r. The wave propagates out of the plane in (a), (b), (c), 
and it propagates in the plane in (d), (e) and (f). Where 
in (a) and (b) we have respectively the plus mode and cross 
mode, in (c) the scalar mode, in (d), (e) and (f) the D, B and 
C mode. 



and the amplitude can be written in terms of the 6 
polarization states as 

V,(t, z) =A+(t- v a . a z)e$ +A*(t- v a . a z)e$ 
+B B (t - v Gs2 z)e$ + C c (t- v Ga2 z)efJ 
+D D (t - v G . a z)e$ + h s (t - v G z)e%. 

(3.3) 

where vo s2 is the group velocity of the massive spin-2 
field and is given by 



V G S 



(3.4) 



The first two polarizations are the same as in the mass- 
less case, inducing tidal deformations on the x-y plane. 
In Fig.l, we illustrate how each GW polarization affects 
test masses arranged on a circle. 

The presence of the ghost mode may seem as a pathol- 
ogy of the theory from a purely quantum-mechanical ap- 
proach. There are several reasons to consider such a 
mode as problematic if we wish to pursuit the particle 
picture interpretation of the metric perturbations. The 
ghost mode can be viewed as either a particle state of 
positive energy and negative probability density, or a pos- 
itive probability density state with a negative energy. In 



the first case, allowing the presence of such a particle 
will quickly induce violation of unitarity. The negative 
energy scenario leads to a theory where there is no mini- 
mum energy and the system thus becomes unstable. The 
vacuum can decay into pairs of ordinary and ghost gravi- 
tons leading to a catastrophic instability. 

One way out of such problems is to impose a very weak 
coupling of the ghost with the rest of the particles in 
the theory, such that the decay rate of the vacuum will 
become comparable to the inverse of the Hubble scale. 
The present vacuum state will then appear to be suffi- 
ciently stable. This is not a viable option in our theory, 
since the ghost state comes in the gravitational sector, 
which is bound to couple to all kinds of matter present 
and it seems physically and mathematically unlikely for 
the ghost graviton to couple differently than the ordinary 
massless graviton does. Another option is to assume that 
this picture does not hold up to arbitrarily high energies 
and that at some cutoff scale M cuto ff the theory gets 
modified appropriately as to ensure a ghost-free behav- 
ior and a stable ground state. This can happen for ex- 
ample if we assume that Lorentz invariance is violated 
at M cuto ff, thereby restricting any potentially harmful 
decay rates [33j . 

However, there is no guaranty that theories of modified 
gravity such as the one investigated here are supposed to 
hold up to arbitrary energies. Such models are plagued 
at the quantum level by the same problems as ordinary 
General Relativity, i.e. they are non-renormalizable. It is 
therefore not necessary for them to be considered as gen- 
uine candidates for a quantum gravity theory and the cor- 
responding ghost particle interpretation becomes rather 
ambiguous. At the purely classical level, the perturbation 
hp, v should be viewed as nothing more than a tensor rep- 
resenting the "stretching" of spacetime away from flat- 
ness. A ghost mode then makes sense as just another way 
of propagating this perturbation of the spacetime geome- 
try, one which carries the opposite sign in the propagator 
than an ordinary massive graviton would. 

Viewed in this way, the presence of the massive ghost 
graviton will induce on an interferometer the same ef- 
fects as an ordinary massive graviton transmitting the 
perturbation, but with the opposite sign in the displace- 
ment. Tidal stretching from a polarized wave on the 
polarization plane will be turned into shrinking and vice- 
versa. This signal will in the end be a superposition of the 
displacements coming from the ordinary massless spin-2 
graviton and the massive ghost. Since these induce two 
competing effects, this will lead to a less pronounced sig- 
nal than the one we would expect if the ghost mode was 
absent, setting in this way less severe constraints on the 
theory. However, the presence of the new modes will also 
affect the total energy density carried by the gravitational 
waves and this may also appear as a candidate signal in 
stochastic backgrounds, as we will see in the following. 
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IV. GRAVITATIONAL WAVES PROPAGATING 
IN A CERTAIN DIRECTION AND THE 
POSSIBLE DETECTOR RESPONSE 

Let us consider now now the possible response of a 
detector revealing GWs coming from a certain direction. 
It is important to stress that the detector output depends 
on the GW amplitude that is determined by a specific 
theoretical model. However, one can study the detector 
response to each GW polarization without speci fying, a 
priori, the theoretical model. Following [19|, I22H25L |39( 
the angular pattern function of a detector to GWs is 
given by 

F A (Cl) = D:e A (ri), (4.1) 

1 



D 



u <K) u - 



here A = +, x , B,C, D, s. The symbol : is contraction 
between tensors. D is the detector tensor representing 
the response of a laser-interferometric detector. It maps 
the metric perturbation in a signal on the detector. The 
vectors u and v are unitary and orthogonal to each other. 
They are directed to each detector arm and form an or- 
thonormal coordinate system with the unit vector w (sec 
Fig. [2j. fl is the vector directed along the GW prop- 
agation. Eq. (|4.1[) holds only when the arm length of 
the detector is smaller and smaller than the GW wave- 
length that we are taking into account. This is relevant 
for dealing with ground-based laser interferometers but 
this condition could not be valid when dealing with space 
interferometers like LISA. 
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Figure 2: The coordinate systems used to calculate the polar- 
ization tensors and the pictorial view of the coordinate trans- 
formation. 



A standard orthonormal coordinate system for the de- 
tector is 

' n= (1,0,0) 
v = (0,1,0) . 
w= (0,0,1) 



On the other hand, the coordinate system for the GW, 
rotated by angles (9, 0), is given by 

u' = (cos 9 cos 0, cos 9 sin 0, — sin 9) 
v' = (— sin 0, cos 0, 0) 
w' = (sin 9 cos 0, sin 9 sin 0, cos 9) 

The rotation with respect to the angle if), around the 
GW-propagating axis, gives the most general choice for 
the coordinate system, that is 

m = u' cos ip + v' sin ip 
n = — v' sin?/; + u' cosip . 



Coordinates (u, v, w) are related to the coordinates 
(m, n, ti) by the rotation angles (0, 9, V>), as in Fig. [5] 
By thevectors rh, n, and f2, the polarization tensors are 



e+ = 



e x = 



—= (m < 

— = (m< 

V2 



> m — n 



) n + n 



e c 



e, = 



— = [rh®n 

V2 



V2 
2 



n 



m 
~2 



m 
~2 



in) 

m) 

• O (g) m 

Cl (g) n 
n 

2 g 



n 
2 



Taking into account Eqs. (|4.1[) . the angular patterns 
for each polarization are 



F+(0,<f,,il>) 



( 1 + cos 2 9) cos 2(j) cos 2ip 



1 

— cos 9 sin 2(f> sin 2tp 



Fx (9, ' 

F B (9,< 
F c (9,< 

F D { 
FJ 



,1>) = 
,V>) 



cos 9) cos 2(f> sin 2if> 



— cos 9 sin 2(f> cos 2ifj , 

sin 9 (cos 9 cos 20 cos ifj — sin 20 sin ip) , 
sin 9 (cos 9 cos 2(f) sin if> + sin 2(f) cos ip) , 

— cos 2</> (6 sin 2 6< + (cos 29 + 3) cos 2ip) , 
32 



V2 



sin 9 cos 20 . 



The angular pattern functions for each polarization are 
plotted in Fig. [3] These results, also if we have considered 
a different model, are consistent, for example, with those 
in [l9l - |2lj . Another step is now to consider the stochas- 
tic background of GWs in order to test the possible de- 
tectability of such further contributions in gravitational 
radiation. 
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Figure 3: Plots along the panel lines from left to right of 
angular pattern functions of a detector for each polarization. 
From left plus mode F+, cross mode F x , B mode Fb, C mode 
Fc, D mode Fb,and scalar mode F s . The angular pattern 
function of the Fb and Fc mode is the same except for a 
rotation. 



V. THE STOCHASTIC BACKGROUND OF 
GRAVITATIONAL WAVES 



The contributions to the gravitational radiation com- 
ing from higher order gravity could be efficiently selected 
if it would be possible to investigate gravitational sources 
in extremely strong field regimes. In such a case, the fur- 
ther polarizations coming from the higher order contribu- 
tions could be, in principle, investigated by the response 
of a single GW detector described above. However, this 
situation seems extremly futuristic at the moment so the 
only realistic approach to investigate these further con- 
tribution seems the cosmological background, in partic- 
ular, the stochastic background of GWs. Such a GW 
background can be roughly divided into two classes of 
phenomena: the background generated by the incoherent 
superposition of gravitational radiation emitted by large 
populations of astrophysical sources (hard to be resolved 
individually (34|), and the primordial GW background 



generated by processes in the early cosmological eras [35| . 
Primordial components of such background are interest- 
ing, since they carry information on the primordial Uni- 
verse and, on the other hand, can give information on 
the gravitational interaction at that epochs [40j, |4l| . The 
physical process of GW production has been analyzed, 
for example, in [361 - 138} but only for the first two stan- 
dard tensorial components of Eq. (|3.2p . that is the GR 
components. Actually the process can be improved con- 
sidering all the components that we have considered here. 
Before starting with the analysis, it has to be emphasized 
that, considering a stochastic background of GWs, it can 
be described and characterized by a dimensionless spec- 
trum (see the definition [H [H [Hi E| ) 



1 dpjw 
Pc din/' 



where 



8nG 



(5.1) 



(5.2) 



is the (actual) critical energy density of the Universe, Ho 
the today observed Hubble expansion rate, and dp sgw is 
the energy density of the part of the gravitational radia- 
tion contained in the frequency range / to / + df. 



p gw = / dfpgwif) 
lo 



(5.3) 



where /5gw is the GWs energy density per unit frequency. 
n gw (/) is related to S h {f) by [HlH 



/4tt 2 



(5.4) 



Note that the above definition is different from that in 
the literature [H, [39j], by a factor of 2, since it is defined 
for each polarization. It is convenient to represent the 
energy density with the form /iq r2 gw (/) by parametriz- 
ing the Hubble constant as Ho = 100 ho kms -1 Mpc -1 . 
Then, the GW stochastic background energy density of 
all modes can be written as 



Q A = 0+ + Q X 4-Q B 4- Op -L Cl D _|_ o s 
gw — J : gw ~ ^ ^gw ~ a 'gw ~ ^ 'gw ' 2 'gw ' 2 'gw 



(5.5) 



we can split ilg W as a part arising from GR 

n GR = n+ +n x , n+ =n x (5.6) 

gw gw gw ' gw gw V / 

a part from higher-order-gravity 

n c = n D (5.7) 



n HOG = +n C +Q D 



gw u "gw 



and a scalar part fig W . 

We are considering now standard units and study only 
the modes which arise from higher order theory. 
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The relic stochastic background of GWs can be de- 
rived by considering only general assumptions and ba- 
sic principles of Quantum Field Theory and GR. The 
quantum fluctuations of the zero-point energy can be 
amplified in the early Universe by the large variations 
of gravity and this mechanism produces GWs. A very 
interesting by-product of GWs is that they can be used 
to probe the evolution of the Universe at early times, 
eve n up to the Planck epoch and the Big Bang singular- 
ity [11 HE H| |Hj|. The mechanism of the GWs is con- 
nected to inflationary scenario 0, EH, which fits well 
the WMAP data and is in particularly good agreement 
with almost exponential inflation and spectral index rs 1 , 

mm. 

A remarkable fact about the inflationary scenario is 
that it contains a natural mechanism which gives rise 
to perturbations for any field. It is important for our 
aims that such a mechanism provides also a distinctive 
spectrum for relic scalar GWs. These perturbations in 
inflationary cosmology arise from the most basic quan- 
tum mechanical effect: the uncertainty principle. In this 
way, the spectrum of relic GWs that we could detect to- 
day is nothing else but the adiabatically-amplified zero- 
point fluctuations [H, H3| • The calculation for a simple 
inflationary model can be performed for the scalar field 
component of eq. (|3.2[) . Let us assume that the early Uni- 
verse is described an inflationary de Sitter phase emerg- 
ing in a radiation dominated phase [HI, [37], SH- The 
conformal metric element is 

ds 2 = a 2 (ri)[-dr? + dl? 2 + h^ v (r),^)dx^dx' / }, (5.8) 



where, for a purely GW the metric perturbation (J3T2J) 
reduces to 



h A e 



(A) 



(5.9) 



where A = +, x , B,C, D, and s. Let us assume a phase 
transition between a de Sitter and a radiation-dominated 
phase [36l |37| , we have: r/i is the inflation-radiation tran- 
sition conformal time and rjo is the value of conformal 
time today. If we express the scale factor in terms of 
comoving time cdt = a(t)dt], we have 



a(t) oc exp(iJrf s t), a(t) oc y/i 



(5.10) 



for the de Sitter and radiation phases respectively. In 
order to solve the horizon and flatness problems, the con- 
dition —7-^-7 > 10 27 has to be satisfied. The relic scalar- 
»Wi) 

tensor GWs are the weak perturbations h^^, x) of the 
metric (|5 .9[) which can be written in the form 

V =e${k)X{ri)e*v{ii (5-11) 

in terms of the conformal time 77 where k is a constant 
wavevector. From ea. Q5.llJ> . the component is 

^T},le,~£)=X(rj)exp(it-~£). (5.12) 



Assuming Y(rf) — a(rf)X(rf), from the Klein-Gordon 
equation in the FRW metric, one gets 



Y 



k\ 2 )Y = 

a 



(5.13) 



where the prime ' denotes derivative with respect to the 
conformal time. The solutions of eq. f|5 . 13J) can be ex- 
pressed in terms of Hankel functions in both the infla- 
tionary and radiation dominated eras, that is: 
For 77 < 771 



x( v ) 



aim) 

a(rf) 



[l + iH ds uj 1 ]exp(-ifc(?7-?7i)), (5.14) 



for rj > r/i 



X ^) = ^M[aexpH%- 
a(ri) 



m)) +0exp(ik(r)-r) 1 ))] , 



(5.15) 

where 10 — ck/a is the angular frequency of the wave 

(which is function of the time being k = \ k | constant), 
a and /3 are time- independent constants which we can 
obtain demanding that both X and dX/drj are contin- 
uous at the boundary 77 = 771 between the inflationary 
and the radiation dominated eras. By this constraint, we 
obtain 



a = 1 + i 



y/HdsHo Hd s H 
u> 2w 2 



p = ^- (5 - 16) 



In eqs. (|5.16p . u) = ck/a(rio) is the angular frequency 
as observed today, H n = c/r/ is the Hubble expansion 
rate as observed today. Such calculations are referred in 
literature as the Bogoliubov coefficient methods [3(1 H3] ■ 
In an inflationary scenario, every classical or macro- 
scopic perturbation is damped out by the inflation, i.e. 
the minimum allowed level of fluctuations is that required 
by the uncertainty principle. The solution (15. 14J) corre- 
sponds to a de Sitter vacuum state. If the period of in- 
flation is long enough, the today observable properties of 
the Universe should be indistinguishable from the prop- 
erties of a Universe started in the de Sitter vacuum state. 
During the radiation dominated phase, the particles are 
described by the eigenmodes that correspond to the co- 
efficients of a, while the antiparticles correspond to the 
coefficients of j3. Therefore, the number of particles that 
have been created at angular frequency uj in the radiation 
phase is given by 



HdsHo 
2w 2 



(5.17) 



Now it is possible to write an expression for the en- 
ergy density of the stochastic scalar-tensor relic gravitons 
background in the frequency interval (uj, uj + doS) for each 
mode as 



dpi 



huj 



2tt 2 c 3 



hH 2 s H 2 duj hH 2 ds H 2 df 



87T 2 C 3 UJ 



8tt 2 c 3 / 



(5.18) 
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where /, as above, is the frequency in standard comoving 
time. eq. (|5.18l) can be rewritten in terms of the today 
and de Sitter value of energy density being 



8irG Pc 
3c 2 



8nGp, 



ds 



3c 2 



Introducing the Planck density ppi a 
trum is given by 



KG 2 



yw 



1 dp. 



yw 



f dp 



yw 



Pds 



Pc d\nf 



df 9 ppi 



anck 



(5.19) 



the spec- 



(5.20) 



At this point, some comments are in order. First of all, 
such a calculation works for a simplified model that does 
not include the matter dominated era. If we also include 
such an era, we would also have to take into account the 
redshift at the equivalence epoch and this results in [38| 



Pds 



9 PPla 



-{i + z eq y 



(5.21) 



for the waves which, at the epoch in which the Universe 
becomes matter dominated, have a frequency higher than 
H eqi the Hubble parameter at equivalence. This situa- 
tion corresponds to frequencies / > (1 + z eq ) x l 2 H§. The 
redshift correction in eq. (|5.21j) is needed since the to- 
day observed Hubble parameter Hq would result differ- 
ent without a matter dominated contribution. At lower 
frequencies, the spectrum is given by [36|, 



(5.22) 



As a further consideration, let us note that the results 
(|5.20|) and (|5.21j) . which are not frequency dependent, 
do not work correctly in all the range of physical fre- 
quencies. Waves that have frequencies less than Hq, the 
energy density is in a sense not well defined, as their 
wavelength becomes larger than the Hubble scale of the 
Universe. In a similar manner, at high frequencies, there 
is a maximal frequency above which the spectrum rapidly 
drops to zero. In the above calculation, the simple as- 
sumption that the phase transition from the inflationary 
to the radiation dominated epoch is instantaneous has 
been made. In the physical Universe, this process occurs 
over some time scale At, being 



fmax — 



a(ti) 1 
a(t ) At' 



(5.23) 



which is the redshifted rate of the transition. In any 
case, Sl£ w drops rapidly. The two cutoffs at low and high 
frequencies for the spectrum guarantee that the total en- 
ergy density of the relic gravitons is finite. These results 
can be quantitatively constrained considering the recent 
WMAP release. Nevertheless, since the spectrum falls 
off oc f~ 2 at low frequencies, this means that today, at 
LIGO- VIRGO and LISA frequencies, one gets for the GR 
part [39] [42] 



< 2 x 10" 



(5.24) 



for the higher-order-gravity part 

^C G (/)^oo<6.7xlO- 9 . 
and for the scalar part 

n™(/)^ioo < 2 -3 x 10- 12 . 



(5.25) 



(5.26) 



It is interesting to calculate the corresponding strain at 
~ 100Hz, where interferometers like VIRGO and LIGO 
reach a maximum in sensitivity [f| 0]- With a minor 
modification we can use the well known equation for the 
characteristic amplitude [39[ for one of the components 
oftheGWs 3 : 



h A {f) ^ 8.93 x 10 



_ 19 flHz 



V / 



V^oo^W/), (5-27) 



and then we obtain for the GR modes 

h GR (100Hz) < 1.3 x 10~ 23 . (5.28) 

while for the higher-order modes 

h HOG {l00Hz) < 7.3 x 10~ 25 . (5.29) 

and for scalar modes 

h s (100Hz) < 2 x 1.410" 26 . (5.30) 

Then, since we expect a sensitivity of the order of 10~ 22 
for the above interferometers at s» 100-ffz, we need to 
gain at least three orders of magnitude. At smaller fre- 
quencies the sensitivity of the VIRGO interferometer is 
of the order of 10 -21 at sa 10Hz and in that case it is for 
the GR modes 

h GR (100Hz) < 1.3 x 10~ 22 . (5.31) 

while for the higher-order modes 

h HOG {W0Hz) < 7.3 x 10~ 24 . (5.32) 

and for scalar modes 

h s (W0Hz) < 1.4 x 10" 25 . (5.33) 



Still, these effects are below the sensitivity threshold to 
be observed. The sensitivity of the LISA interferometer 
will be of the order of 10" 22 at 10~ 3 iiz (see [|) and 
in that case it is 



h GR (100Hz) < 1.3 x 10 



-18 



(5.34) 



3 The difference between our result and eq. (19) in Ref. [39l | is due 
to the fact that the latter did their calculation assuming the two 
polarization modes of GR while we handle each mode separately, 
hence the difference. 
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while for the higher-order modes 

h HOG (100Hz) < 7.3 x 10~ 20 . (5.35) 
and for scalar modes 

h a (W0Hz) < 1.4 x 10~ 21 . (5.36) 

This means that a stochastic background of relic GWs 
could be, in principle, detected by the LISA interferom- 
eter, including the additional modes. 

VI. CONCLUSIONS 

Our analysis covers extended gravity models with a 
generic class of Lagrangian density with higher order and 
terms of the form f(R,P,Q), where P = R a bR ab and 
Q = Rabcd R abcd . We have linearized the field equations 
for this class of theories around a Minkowski background 
and found that, besides a massless spin-2 field (the gravi- 
ton), the theory contains also spin-0 and spin-2 massive 
modes with the latter being, in general, ghosts. Then, 
we have investigated the detectability of additional po- 
larization modes of a stochastic GW with ground-based 
laser-interferometric detectors and space-interferometers. 
Such polarization modes, in general, appear in the ex- 
tended theories of gravitation and can be utilized to con- 
strain the theories beyond GR in a model-independent 
way. 

However, a point has to be discussed in detail. If the 
interferometer is directionally sensitive and we also know 
the orientation of the source (and of course if the source is 
coherent) the situation is straightforward. In this case, 
the massive mode coming from the simplest extension, 
/(i?)-gravity, would induce longitudinal displacements 
along the direction of propagation which should be de- 
tectable and only the amplitude due to the scalar mode 
would be the true, detectable, "new" signal [27]. But 
even in this case, we could have a second scalar mode 
inducing a similar effect, coming from the massive ghost, 
although with a minus sign. So in this case, one has de- 
viations from the prediction of /(i?)-gravity, even if only 
the massive modes are considered as new signal. 

On the other hand, in the case of the stochastic back- 
ground, there is no coherent source and no directional 



detection of the gravitational radiation. What the inter- 
ferometer picks is just an averaged signal coming from 
the contributions of all possible modes from (uncorre- 
cted) sources all over the celestial sphere. Since we ex- 
pect the background to be isotropic, the signal will be the 
same regardless of the orientation of the interferometer, 
no matter how or on which plane it is rotated, it would 
always record the characteristic amplitude h c . So there 
is intrinsically no way to disentangle any of the mode in 
the background, being h c related to the total energy den- 
sity of the gravitational radiation, which depends on the 
number of modes available. Every mode, essentially, con- 
tributes in the same manner, at least in the limit where 
the mass for the massive and ghost modes are very small 
(as they should be). So, it should be the number of the 
modes available that makes the difference, not their ori- 
gin. 

Again, even if this does not hold, one should still get 
into consideration at least the massive ghost mode to 
get a constraint. This is the why we have considered 
only hen, htfOG and h s in the above cross-correlation 
analysis without giving further fine details coming from 
polarization. For the situation considered here, we find 
that the massive modes are certainly of interest for direct 
attempts at detection with the LISA experiment. It is, 
in principle, possible that massive GW modes could be 
produced in more significant quantities in cosmological 
or early astrophysical processes in alternative theories of 
gravity, being this possibility still unexplored. This situ- 
ation should be kept in mind when looking for a signature 
distinguishing these theories from GR, and seems to de- 
serve further investigation. 
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